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Abstract 

We consider the simple hypothesis of letting quantum systems have an 
inherent random nature. Using well-known stochastic methods we thus 
derive a stochastic evolution operator which let us define a stochastic 
density operator whose expectation value under certain conditions satisfies 
a Lindblad equation. As natural consequences of the former assumption 
decoherence and spontaneous emission processes are obtained under the 
same conceptual scheme. A temptative solution for the preferred basis 
problem is suggested. All this is illustrated with a comprehensive study 
of a two-level quantum system evolution. 



1 Introduction 

Stochastic methods arc lately being profusely used in Quantum Mechanics in 
general and in Quantum Optics in particular to better study and analyze the 
physical processes taking place during the interaction between the system and 
the measuring apparatus. The conceptual framework is very clear, consisting 
in an open quantum system which is being continuously monitorized by an 
appropiate device. This situation is very well described using stochastic un- 
ravellings for the Lindblad equation of the system density operator, the most 
general expression of which for the diffusive case has been recently given in Q . 
Needless to say, the presence of an observer is compulsory for orthodox quantum 
principles to be applied. 



On the other hand, in the early years of Quantum Mechanics, the role of the 
observer was crucial to both justify the indeterministic nature of the results of a 
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measurement upon quantum systems g and to understand the lack of quantum 
interference in a, say, double-slit experiment. Nowadays it is understood that 
this lack of interference is not strictly due to the presence of the observer, but 
to the possibility of knowing which path the particle followed ||. Keeping 
orthodox principles Q as far as possible, the role of the observer is then left as 
the only justification to set the indeterminism in the measurement process of 
quantum systems. 

In this article, we take under consideration the possibility of reducing once 
more the role of the observer and study the consequences of letting quantum sys- 
tems have an intrinsic stochastic nature. In the same fashion, this hypothesis has 
already given rise to different stochastic evolution models (cf. ||] and references 
therein to previous works), the main purpose of which was to reproduce faith- 
fully all nonrelativistic quantum predictions with the sole exception of quantum 
superpositions among macroscopically distinguishable states. By concentrating 
upon a very simple quantum system, namely a two-level system, we will show 
how this stochastic assumption allows us to make predictions only restricted 
to relativistic extensions of Quantum Mechanics (QED) or only achievable by 
resorting to the theory of open quantum systems. We will restrict ourselves to 
the nonrelativistic domain, i.e. we will keep on using state vectors to describe 
physical systems, something which is impossible to do in an attempt to merge 
relativistic principles with Quantum Theory. 

2 Motivation and General Framework 

The main physical hypothesis we address here is the possibility of endowing the 
evolution operator of a quantum system with stochastic nature. To be concrete 
we will deal with a two-level quantum system in particular. This assumption 
has already been done in another places (cf. e.g. |6) for a general overview) 
and can be alternatively motivated in two ways. On one hand any system, 
whether quantum or classical, is continuously subjected to the external influence 
of its environment, however small the said influence may be. This standpoint 
is assumed in the program of decoherence, in fact it is its conceptual starting 
point Such external influences are, needless to say, uncontrollable, so we 
may encode them using stochastic methods. This way of proceeding is similar 
to for instance the modelling of the evolution of a Lorentz particle Qj. In van 
Kampen's terminology it corresponds to what he calls external noise which is 
due essentially to the huge amount of uncontrollable external factors that affect 
the system evolution. 

Alternatively, we may adopt the complementary point of view, stating that 
the stochasticity arises just as internal noise (cf. again 0). This second stand- 
point is subtler and, in our opinion, much less intuitive but hopefully more 
accurate. Let us use a canonical example (which will be extensively studied 
below), a two- level quantum system in electromagnetic vacuum. Within the 
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non-relativistic quantum-mechanical standpoint the evolution of this system 
is completely determined by its hamiltonian Hq and since the latter is time- 
independent, the system is stationary, i.e. the energy does not change. Within 
QED framework, the description varies to include the effect of vacuum fluctu- 
ations. These fluctuations, quantum in nature, drive the system to its ground 
energy level, thus producing spontaneous decay. This situation can also be un- 
derstood letting the system be modelled by stochastic evolution, rooting the 
stochasticity in this quantum vacuum. Since these fluctuations do not occur in 
real ordinary space, but in Hilbert space we use stochastic evolution operators 
to describe them. The term internal is in this situation a bit dubious, since 
the electromagnetic vacuum can hardly be thought as internal to the two-level 
system. We should understand then internal as making reference to the essence 
itself of the system. As a matter of fact this vacuum is only detectable through 
its effects on a quantum system, so in certain sense the evolution of the system 
should also contain such effects. 

From a mathematical point of view the question of the origin of the stochas- 
ticity is secondary and eveything reduces to find the form of such a stochastic 
evolution operator. So we proceed by substituting U(t) U s t(t) and then in- 
vestigating the form of U st {t). To do this we resort to a general decomposition 
theorem of real random variables @] the generalization to (bounded) operator- 
valued random variables of which we take for granted. 

Theor. 2.1. Let X be an operator-valued random variable acting upon a N- 
dimensional Hilbert space Then there exist N x N operator-valued processes 
Vk{t) such that 



where E denotes the expectation value with respect to the probability measure 
and Wt are n complex Wiener processes []. 

Expressing the latter integral as a function of the upper interval limit (T 
R + ) we may write as the more general form for a stochastic evolution operator 



where in general the Wiener processes will have the covariance matrix given 




(2.1) 




(2.2) 



by 



E[W t k W t k '*} = Ofcfc/t. 



(2.3) 



iCf. e.g. 



for a general reference on stochastic methods. 
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The connection with ordinary Quantum Mechanics is made by stating 



p QM (t)=Ep st (t), (2.4) 

where p s t{t) is the system density operator induced by the previous stochas- 
tic evolution operator, namely 



Pst (t) = U st (t)p(0)U f st (t). (2.5) 

Henceforward, we adopt the same assumptions as in the orthodox formalism, 
in particular we force trace-preserving evolution. This immediately imposes 
conditions on the absolute value of ¥,U s t(t), as the following proposition shows: 



Prop. 2.1. If pQiait) is defined as and trpQM{t) — trpQM(0), then 

1/2 



EU st (t) = exp(-iHt) 7- a nm / uj n (s)u„(s)d 

\ nm=l Jo 



s\ , (2.6) 



where H is a selfadjoint operator. 

Proof. Cf. Appendix A. □ 

A few comments should be included. Firstly the selfadjoint operator H is to 
be identified with the hamiltonian of the system. Secondly note, as expected, 
that if v(t) = 0, then we recover the usual quantum-mechanical formalism, thus 
reinforcing the idea that new effects are contained only in the stochastic part 
of the evolution operator, leaving the hamiltonian partial evolution unchanged. 
Thirdly the square root imposes restriction upon J2 nm a nmvJ n (t)v n (t) since 

I- anm vi(s)v n (s)ds > (2.7) 

nm=l ^° 

must hold for all t > 0. Finally U s t{t) is not an unitary operator, unitarityQ 
is only obtained in stochastic average, i.e. 



E[Uj t (t)U st (t)]=I. (2.8) 

This last feature is the crucial difference with other stochastic evolution 
schemes showing the same philosophy (specially those assumed in |l(J and ) . 
There, as in the scheme we propose here, it is assumed a basic stochastic evolu- 
tion described by two operators, namely the Hamiltonian H (for the determin- 
istic part) and an ad hoc operator L which gives rise to the Lindblad operator 

2 We are aware that the following condition only refers to partial isometry and not unitarity 
itself. We, for the time being, do not care about mathematical refinements. 
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[ fL2| after the whole calculation is carried on. Nevertheless the stochastic nature 
appears only as a modification of the evolution operator generator. For instance 
example 30. f of jlO] proposes the ansatz 

T(t) \p(0)] = E[e? WtL p(0)e- iWtL ] , (2.9) 
thus obtaining a Lindblad-type generator^ 



9[X] = -{[LX,L\ + [L,XL]}. (2.10) 

This is straightforwardly generalized to more than one Lindblad operator Lk ■ 
In a simpler and more intuitive form much the same it is arrived at in (nfl . They 
only confer stochastic nature to the evolution operator by adding random parts 
to its generator, thus obtaining some further restrictions on the final evolution, 
such as only selfadjoint Lindblad operators and consequently not achieving the 
desirable most general situation (any Lindblad operator). Parthasarathy has 
however provided a scheme to obtain general Lindblad evolution, namely he has 
proposed the evolution 



T t [X] = E[Ui N ®XU N ®], (2.11) 

where N(i) is Poisson process with intensity A and L = \/~\U. Needless to say, 
though he arrives at a Lindblad generator | {[L-,L^] + [L, the physical 

interpretation of such a stochastic evolution operator is rather evasive. So we 
find that to keep the intuitive approach of Adler, we must sacrifice some gen- 
erality whereas on the other hand to arrive at such generality as Parthasarathy 
we lose physical intuition. In order to conjugate both characteristics, we have 
assumed a less restrictive position. 



However this generality does not drive us directly to a Lindblad evolution 
as we shall see. Some restrictions should be made on the added random parts 
Vk(t). These restrictions possess, as the whole scheme, a rather intuitive physical 
interpretation. To see how the Lindblad evolution appears we focus on the 



differential version of (2.4). The starting point is the following 
Prop. 2.2. If p(t) is defined as ( \2.$ ), then it satisfies the differential equation 

p(t) = Lt[p(t)]+L t [p(0)}, (2.12) 

where for any X 

3 The lack of a hamiltonian part can be viewed as a sort of hamiltonian-interaction picture 
election. 
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L t [X] = -i[H,X} + 

^ N 2 

+ g E a nro ([£„(t)X,4(t)] + [^(f),x4(t)]) 

(2.13a) 

Lt[X] = - V a nm L t [/ u n ( a )*t4(«)d8] (2.13b) 
i •/ o 

and by construction we have defined 

e n (t)=v n (t)E[U st ]- 1 (t). (2.14) 
Proof. This result is straighforwardly obtained by deriving with respect to time 



the expression of pqm(^) (2.4) in which we use expression fl2.q ) for U s t(t) and 



identifying the previous definitions. □ 



The equation ( 2.12 ) is almost a Lindblad differential equation. There are 
only two differences: 

1. There is an extra term L t [p(0)] which spoils the Markovianity. This on 
the other hand is also obtained in the more general quantum evolution 
of a subsystem in the orthodox formalism (cf. 1 1 3 or |H). The physical 



conditions to be met to guarantee Markovianity must be found. 

2. The would-be Lindblad operators £ n (t) are time-dependent. Though time- 
dependent Lindblad operators can be found in the literature (cf. |L4| and 
references therein), this is not the usual case upon which we will focus. 

Thus we must find the conditions under which on one hand L t [p(0)] = holds 
and on the other hand £k(t) ■>*» £u- ln the next section we will comprehensively 
analyzed a two-level system to gain more physical insight into the question 
instead of trying to obtain the general conditions. 

However notice that a first partial general result can be obtained as to the 
energy preserving or not preserving nature of the fluctuation part v(t) once the 
second term L t [p(0)] is shown to be negligible: 

Prop. 2.3. Both [v(i),v^ (t)] — and [H,v(t)} = are independent sufficient 
conditions on their own to get energy conservation, i.e. 



d -^ = ltrHp(t)=0. (2.15) 
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Proof. Its elementary using ( |2.12 ) with i t [p(0)] = and the cyclic property of 
the trace. □ 



This result may be useful to check whether we have a decaying or pumping 
interaction or on the contrary we just have a decohering but energy-preserving 
evolution (see below). 

Before working out the announced two-level quantum system example, we 
will briefly state where the main differences between previous stochastic quan- 
tum evolution models and our scheme are. Let us consider the Quantum State 
Diffusion model (cf. JTHJ ) . This model corresponds to stochastic unravellings of 
the Lindblad equation under the hypothesis of state- vector normalization, which 
is implemented in the following way. If the stochastic state-vector satisfies the 
QSD equation 

d\ip st ) = \v)dt + \f)dt, (2.16) 
then the normalization preservation is achieved by forcing 

(^t|/> = 0. (2.17) 
As \ij} s t) must stochastically unravel the Lindblad equation, then 

pit) =n\i> st {t))^ st {t)\). (2.i8) 

The QSD unravelling for the Lindblad equation with one Lindblad operator 

is 



d\i> 8t ) = -iH\i> st )dt + (V)L - h^L - ~(tf)(L)\ \ip st )dt + (L- (L))\rp st )d^ 

(2.19) 

where L denote the Lindblad operator entering the corresponding Lindblad 
equation and £ denotes white noise (cf. [^5|). Note that this equation is non- 
local and non-linear. The difference with the scheme proposed here is double. 
Firstly we do not assume the Lindblad equation as a starting hypothesis, we 
arrive at it only after finding some physical conditions (cf. below) . Secondly we 
do not impose normalization of \ip s t), we only impose the condition 



n\(At(tMst(t))\ 2 } = ^ 



(2.20) 



i.e. normalization in stochastic average. To better appreciate these differ- 
ences, we will explicitly write down the evolution equation for \ipst) derived from 
the evolution operator (2.2): 
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d\At) = (-iH - -£\t)l{t)W st )dt 
+ v(t)\i>(0))dW t . 



(2.21) 



Notice that this equation is linear and, under appropiate elections of v(t), 
local. Notice also the explicit appearance of IV'(O)), in consonance with the most 
general evolution equation for an open quantum system JL3| . But QSD is not 
the unique diffusive stochastic unravelling for the Lindblad equation (cf. Q). 
As a matter of fact it is included as a particular case of the stochastic evolution 
equation given by Diosi and Wiseman. Again, in their approach the Lindblad 
equation is assumed from the very beginning, which settles the first difference 
with the formalism proposed here. Secondly, as a consequence of the assumed 
conceptual framework (a quantum system being continuously monitorized) , they 
impose the condition of purity conservation for the stochastic evolution, which 
gives rise to the main analytical difference, since in our approach U s t(t)'U s t(t) ^ 
/, as we hav e indicated above. Unitarity is only satisfied in stochastic average 



Differences can also be found with the work of Ghirardi, Pearle and Rimini 
[|]. There the Lindblad equation is not assumed from the beginning, but it 
is arrived at after imposing some other conditions. Among them, the most 
important one is, in our opinion, the election of the quantity to be containing the 
physical information of the quantum system. Their starting stochastic evolution 
equation is linear and (if appropiate elections are made) local, just as in our 
approach: 



but they claim that, since this evolution does not preserve normalization, 
the physical information cannot be contained in So they define another 
stochastic vector \(f>) by 




(2.23) 



with probability function given by the probability function corresponding to 
\ip) times their squared norm ||^>|| , i.e. 



where q{<f)) denotes the probability associated to \<f> s t) an d p(i/j) denotes 
the probability associated to \ip s t)- This assumption leads immediately to the 
nonlinear stochastic differential equation 



(cf. eq. fl2.8|)). 



d\il>) = C\*p)dt + A\ip) -dB, 



(2.22) 



g(*) = IMI 2 p(V0, 



(2.24) 
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d\4>) 



iff- J(At-H^).A + J(A-R^ 



)dt + (A - 



(2.25) 



where 7 denotes the non-null elements of the diagonal covariance matrix of 
the set of Wiener processes B and = ^(</>|A + AT|^>). Now the difference 
is rooted in the assumption that in our approach nonnormalization does not 
suppose a major nuisance, since we claim that the physical information of the 
system is contained in 



p(t)=E[|Vu><lM], (2-26) 
so that the probability of finding a system in the state P a = \o)(cr\ will be 



trpP CT =E[|ftM<r)| 2 ]. (2.27) 

And it is the physical quantity ( |2.26| ) which is to be normalized. Thus we 
do not need to resort to some other normalized process \<p). 

Finally, though being very close in spirit, some differences can also be re- 
marked with the work of Gisin 0. Here the objective is to classify all the pure 
state valued stochastic differential equations in C 2 such that the corresponding 
density matrix follows a quantum dynamical semigroup evolution. The first 
formal difference stems form the fact that he describes the stochastic evolution 
in terms of the "height" and azimutal angle of the corresponding point in the 
Bloch Sphere, thus preventing its analysis from being generalized to higher-level 
systems in a straightforward fashion. Moreover he only focus on pure states in 
contrast to our density matrix formalism, which also embraces mixture states. 
Last of all, again as in previous commented models, the evolution equation for 
the density operator is assumed from the beginning, though the restriction of 
complete positivity is not guaranteed in his analysis. 



3 Non-perturbed Two-level Systems 

We begin by considering non-perturbed systems, i.e. with hamiltonian of the 
form 



flb-(t I) (3-28) 

The reason to focus only on non-perturbed hamiltonians will be clear in the 
next section. As usual (cf. e.g. |]l6)) we may choose E + + E- = 0, thus we 
write 
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H = fiuoQ (3.29) 

As a canonical example we may consider an atom in electromagnetic vacuum. 
The random part of the stochastic evolution operator may be written in the 
energy eigenvector basis as 



J*Xn{a)dWi 0\ (0 J*\ 12 (s)dW? 

0/ v° 

0\ /0 

X\2i(s)dW? o) + [o !l\^[s)dWt 

This is the most general form. To clearly understand the physical meaning 
of each of the previous four terms we will start by considering them one by one. 



(3.30) 



3.1 Single Couplings: Decohering, Decaying and Pump- 
ing Factors 

Thus let us study the evolution induced by a stochastic evolution operator of 
the form 



U st {t) 



Hint 



v' (s)v(s)ds 



1/2 



v(s)dW s , 



(3.31) 



where the stochastic expectation value has already been written to preserve 
the trace, v(t) is one of the previous Xij(t)Eij's {{Eij} being the canonical basis 
of A^2(C)) and Wt is the complex standard Wiener process. As a first result we 
elementarily arrive at L t [p(0)] = for any election of v(t), so Markovianity is 
straightforwardly achieved. To obtain time-independent Lindblad operators, we 
insert sucessively each v(t) in the definition of £(t) (eq. (2.14)). Forcing time- 
independency sets an ordinary differential equation for |Aij(i)| whose solution 
leaves as the only option for Ay (t) the expression 



It 



=F iujQt 



(3.32) 



The Lindblad operators produced in this way are of the form 



t = 7 £y (3.33) 

The physical interpretation is rather intuitive. If for instance we choose 
v(t) = Xn(t)Eii, then we obtain an energy-preserving decohering evolution^ , 

4 As usual pij(t) denote the ijth entry of the density matrix p(t). 
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i.e. 



_ = --[F ,p(t)] + - 











(3.34) 



This means that W^i) represents a vacuum-system interaction which does 
not change the energy of the system, but that introduces decohering factors. 
Notice that exactly the same equation (3.34) would have been obtained had 
we chosen v(t) = \22(t)E22- The difference stems from the fact that whereas 
W l {i) represents an interaction through coupling to the excited energy level, 
W A (t) represents a similar interaction but coupled to the ground energy level. 
Graphically the situation can be represented as in Fig. [y. 



Excited-Level-Coopled Decohering Factor 



■ Ground-Level-Cotipled Decohering Factor 



Figure 1: Representation of energy-preserving decohering factors in the stochas- 
tic evolution of two-level quantum system. 



This same interpretation can be carried over to the other two terms. For 
instance, for v(t) — A2i(£)£ , 2i the final evolution equation is 



dp(t) 
dt 



[H , p{t)] 



-2pm(t) -Pu(t) 
-Pai(*) 2/Jaa(*) 



(3.35) 



Figure 2: Spontaneous Decay of a Two-level Atom. 



Now besides the off-diagonal decohering terms a new effect appears: the 
system spontaneously decays into the ground state. This behaviour is indeed 
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Stochastic Term 


A«(t) 


Lindblad Operator 


Au(i) 


V Oy 




A u (t)= 7 V2 e xp(-f -iwot) 


7 


V Oy 




Aia(i) 


fo 1> 
vO oy 




A12W = 7 1 /2 e xp(-f + two*) 


7 


fo 1\ 
VO oy 




A 21 (i) 


A) <a 

V l Oy 




A 21 (t)= 7 1/2 exp(-f -iw i) 


7 


fo o\ 
v i oy 




A 22 (i) 


A) o^ 

V ly 




A22W =7 1/2 exp(-f +iw t) 


7 


fo o\ 
V ly 





Table 1: Lindblad Operators and Stochastic Parts for single vacuum-system 
couplings. 

exactly what is obtained by orthodox means (cf. e.g. Q and Fig. ||). The 
reverse behaviour, as expected, is achieved if we choose v(t) = Ai 2 (t)£ ! i 2 . In 
these two cases then the vacuum-system interaction is understood as an energy- 
non-preserving (decaying or pumping) interaction, depending upon whether the 
coupling takes place through the excited or the ground state. Pictorially we 
may represent these effects as in Fig. [| The whole results are collected in table 
[l] for each of the previous stochastic parts. 



1 



Excited-Level-Coupled Decaying Factor ■ Ground-State-Coupled Pumping Factor 

Figure 3: Representation of energy-nonpreserving decohering factors in the 
stochastic evolution of two-level quantum system. 



Notice that Markovianity is obtained with no further restrictions upon the 
form of the stochastic parts Vk(t)'s. The final form of their entries is calculated 
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from the time-independence condition for the Lindblad operators. On the con- 
trary, some physical conditions are to be met when two Wiener processes are 
present. 

3.2 Double couplings 

The next natural step is to include one more Wiener process in the stochastic 
evolution, so we must investigate the options 



U st (t) = e~ lHot 



2 ~t 

V] ay / vj(s)vi(s)ds 



1/2 



+ / Vl (a)dW}+ / v 2 {s)dW ! 



(3.36) 



where Vj(t) will alternately be each Xij(t)Eij. The calculations proceed ex- 
actly in the same spirit as before, with the exception that now the covariancc 
matrix EfW^WT™*] = a nm t must be taken into account. The results are con- 
tained in table 0. 



Some comments are to be made. First the impossibilty of having a double 
coupling to the same level is clearly shown (rows 2 and 5), which is quite nat- 
ural from a physical point of view. Perhaps the most appealing feature is the 
necessary uncorrelation between the two Wiener processes expressed through 
the condition ay = 0. There's an immediate interpretation of this condition 
in the basis of Markovian evolution. Physically Markovianity is characterized 
by the fact that the evolution of a system from one initial time to a final time 
is exactly the same as the evolution from that initial time to an intermediate 
time and then from the latter to the final instant, whatever that intermediate 
time is. This strongly suggests the idea that the evolution process is made up of 
very tiny and identical contributions of infinitesimal duration, so that the final 
evolution is just the whole contribution of each infinitesimal part. The common 
view states that such a system lacks of memory, i.e. the state to which it evolves 
only depends on its present state and never in the previous ones. These two 
standpoints are complementary. In our case if the Wiener processes partially 
driving the system evolution were correlated, then this evolution partition into 
tiny equal parts would be impossible, thus neglecting the possibility of Marko- 
vianity. The final form of the matrix entries Ay(t)'s are obtained requiring 
time-independence as before. 



3.3 Triple and Quadruple Couplings 

The considerations taken up before lead us immediately to neglect any possibil- 
ity of establishing triple or quadruple couplings, since each energy level can only 
interact with the electromagnetic vacuum through one single Wiener process. 
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Stochastic Terms 


&ij > Xij (t) 


Lindblad Operators 




v °y 




^0 A 12 (^ 

v° y 




a 12 = 

An(t) - VtT cx P (-^i - «•*>*) 

A 12 W — y 72 cajj ^ 2 t T iLUQUj 


7i 


(i S 


72 








(Xu(t) 0\ 

K °y 




( o o> 

\X 21 {t) Oj 






It's impossible to have both a 
preserving and a nonpreserving 
factor coupled to the same level. 




(Xu(t) 0\ 

v °y 




(0 >| 
\0 X 22 {t) y 




ftl4 = 

An (t) = VtI exp (-^^ - two*) 

^221^ — y74 exp ^ 2 + 'Wo^ 




(J S 


74 








V u u ) 




( o o^ 

\A 2 l(l) UJ 






It's impossible to obtain a 
Markovian evolution 
for any initial p(0). 








fo o A 

\0 X 22 (t)j 






It's impossible to have both a 
preserving and a nonpreserving 
factor coupled to the same level. 




( 0\ 




fo o A 

\0 X 22 {t)j 




a 3 4 = 

A 2 i(i) = 7a 1/2 exp _ j^) 
X 22 (t) =7 4 1/2 cxp(-^+2^) 


73 




74 


(2 5 





Tabic 2: Lindblad Operators and Stochastic Parts for double vacuum-system 
couplings. 
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Only in case of more complex systems we may combine in a variety of different 
ways triple, quadruple and higher-order couplings. 



4 Perturbed Systems 

The generalization of the previous formalism to perturbed, i.e. interacting sys- 
tems poses relevant questions as to the adopted conceptual scheme. Naively we 
may try to investigate the evolution produced by an operator of the form 



U st (t) = Te- i ti H{s)ds 



v{s)dW s 



v\s)v{s)ds 



1/2 



(4.37) 



where H (t) may be nondiagonal, as e.g. that expressing the interaction with 
an electromagnetic wave^] 



H(t) 



hup £e~ 2UJt 
£*e tult -hiu Q 



(4.38) 



Hopefully this as before should drives us to a differential equation for the 
density operator with the form 



^ = -i[H(t),p(t)] + \ {[£p(t),H] + [i,p(0]} . (4.39) 

In particular one should be able to reproduce Rabi oscillations with decay 
(cf. e.g. 17 ), but this is impossible as the following proposition states: 

Prop. 4.1. Let X G A^C) be of the form X — vpv^ where p S A^C) is a 
density matrix and v £ A^2(C) is arbitrary. If in an arbitrary common basis, 
H and £ are expressed as 

H 

then one of the following necessary conditions to have 

L[X] ee -i[H,X] + \ {\IXJ\ + [t,X£i]} = 

H selfadjoint, I arbitrary for all p must hold 

5 The rotating wave approximation has already been assumed. 



All 


h 12 \ 


i = 


[hi 


£l2\ 


U21 


h 2 2 I 
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1. Ifv=(h then 

£ 12 t n = i2h 12 £ ai = (4.40a) 



2. Ifv = 





Pi fa 



then 



£ 21 t 22 = i2h 21 



"12 



= 



(4.40b) 



3. Ifv 



fa 

fa 



fa 
fa 



Pk 's not corresponding to the previous cases, then 



aj.An 

«lAi2 



a 2 A 2 i 
a 2 A 22 



a3A 3 i 
a 3 A 32 



■ CK4A42 



(4.40c) 



where the cti 's and the Ay 's are complicated expressions concerning the 



elements of H and £ (cf. Appendix^). 
Proof. Cf. Appendix 



□ 



To get a decay process we only need case 2- eq. (1.40b)- since the Lindblad 
operator now re ads £ = which is only possible if v is of type ( ^ $ 2 ) 

which following ( 4.40b ) makes the Lindblad components time-dependent and 
the Hamiltonian diagonal. 



Far from invalidating the whole formalism this, in our opinion, suggests the 
correct interpretation to be assumed in the introduction of the stochastic nature 
of the evolution operator: the former scheme must only be applied to closed 
systems. The question then is why is it possible to apply this formalism only 
to traditionally closed systems? The answer can be inspired in the canonical 
example used so far. In an orthodox description, the presence of just a single 
photon changes the whole picture, since now there's a real interaction between 
the electromagnetic field and the atom. In the former case, with no photon, 
i.e. with electromagnetic vacuum, the notion of interaction does not appear 
in a natural way, since there's nothing the atom interacts with. But this does 
not mean that an open (interacting) system may not show stochastic behaviour, 
what this means is that to obtain the correct stochastic evolution operator for an 
open system, we must apply the previous formalism to the system+environment 
compound and then trace out the latter's degrees of freedom. In the general case 
this is nearly impossible since it supposes to be able to control all environmental 
degrees of freedom. To illustrate this idea and check whether the qualitative 
behaviour described in the previous section is still valid, we will work out a 
simple and manageable example, namely a two-level atom described by a Jaynes- 
Cummings Hamiltonian 



16 



Hjp = ojqS z + uja'a + e(a* S- + aS + ). 



(4.41) 



This represents an atom interacting only with one single mode of the electro- 
magnetic field. We will consider only a multiple energy-preserving interaction, 
i.e. for the whole atom+field system we consider a stochastic evolution operator 
of the formal 



U st (t) = EU st (t) + J2^+i°Wt +1 (t)\ut +1 )(ut l+1 \ + 



71=0 



A oW (t)\u )(u |, 



(4.42) 



where A,^(t) = v^T^exp 
form given by Prop. 2.1. 



and EU st (t) has the corresponding 



The previous formalism yields the following solution for the density matrix 
entries: 



„±± 



(t) = exp 



rnm 



f mi 



(0) 



otherwise, 



(0) if n 7^ m or differen^igJui) 
(4.43b) 



which is the similar decohering behaviour previously shown by the single atom. 
But we are interested in the atom itself and not in the whole system, so we must 
trace out the field degrees of freedom, i.e. 



p A {t)=tv FPAF {t). (4.44) 

To be concrete we will focus on two particular situations. First we will 
consider an initially correlated or uncorrelated excited atom and n photons 
with energy hui (n ^ 0). Later we will study the case of no photon (n = 0). 
Under these hypotheses the initial atom+field density matrix in the |±, m) basis 
for the first case is given by 



p{-, X;pq) = 5. + 5x+8pn5qn, 



(4.45) 



where by •, x we denote any of the ± signs. The partial trace (4.44) drives 
us through a tedious though elementary calculation to the result 



5 We denote xf oWf(t) 



f*X±(s)dW±(s). 
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PA 



(pk+lk+1 cos2 k+i - (Pfc+ifc+i - Pk+ik+i) cos6» fe+ i sm9 k+1 + p k+lk+1 sin 2 6> fc+ i) 

k=0 

oo 

^ (Pfe+u cos Ok+i sin 6> fe + /ffe+ lfe cos 8 k+ i cos fc - p k + lk sin 0fc + i sin 6» fe - 



fc=i 



-){-| 



Pfe+ifc+i sm f fc+ i cos pfc j + p w cos fi - p 10 sin tfi 

OO 

^ (PS. ! cos 9 k+ i sin 6> fc - p kk+1 sin 6» fc+ i sin fc + p^ +1 cos 6» fe+ i cos 6 k - 

I-X+I + 



fc=l 



P fc fc + i sm y fc+ i cos y fc J + p 01 cos Wi - p 01 sm f i 



(Pk+ik+i sin2 + (Pfe+ife+i + Pfe-fife+i) cos6> fe+1 sin(9 fc+1 + p k+lk+1 cos 2 6> fe+1 ) 
)(-|, 



k=l 



Poo 



where p'* q are the density matrix componen ts in the \u )} basis 



Notice that all these matrix density entries in (4.46) are time-dependent, the 
time dependence being given by ( f4.43a|) and fl4.43b| ). Indeed as a result of 
(4.43a) and (4.43b), and the initial conditions ( 4.45| ) the final atom state after 
a time t 3> T = max{(7i + 7j)~ 1 } has elapsed is a 2 x 2 matrix given in the 
usual |±) basis by 



-)(" 



(4.46) 



p^ » t ) - ( c - 49 -r in4e " +i w - ( 4 - 4? ) 

from where we can explicitly follow the decoherence suffered by the atom, 
and the probabilities of remaining and decaying in the excited and ground states 
respectively exactly coincide with the quantum ones. As initially chosen (in the 
election of the stochastic part), the energy of atom+field is conserved, so we 
only obtain decoherence in the compound system. However the energy of the 

system may vary, as its reduced density matrix shows (pa( ) ^ 0), this energy 

variation being due to the interaction between the field and the atom. 

The second situation in which there's no photon present (n = 0) is included 
to check consistency with the case of a single atom -see section |[ Now a 
quantitative difference should be expected, since we are only taking into account 
one single mode. However the qualitative behaviour, i.e. the spontaneous decay 
process should also be obtained. The same density matrix form for large t is as 
a matter of fact obtained in a similar calculation 
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^^» T )-( cos49i r in49i 2 cos^ t ; si n^ 1 )- ( 4 - 48 ) 

As expected we do not obtain the same quantitative behaviour as in Q3.35| ). 
However a qualitative spontaneous emission is obt ained , since starting from a 
system completely in the excited state (cf. eq. ( 4.45| )) we have found that 
there's nonnull probability of finding it in the ground level even though there's 
no photon present in the field. 



5 Vacuum Interactions vs. Stochastic Nature 

So far we have kept ourselves under the conceptual orthodoxy derived from the 
Quatnum Field Theory formalism in general, and the Quantum Electrodynam- 
ics in particular. But a natural question can be raised as to the conceptual 
framework which may be attached to the previous section scheme. Is it possible 
to found the motivation of use of stochastic methods on some other vacuum- 
independent notions? This is, we are aware, a delicate question, since it reaches 
the very interpretation of Quantum Mechanics itself. But nothing is further from 
our intention than providing a brand new interpretation of Quantum Mechan- 
ics. We only restrict ourselves to consider the possibility of merging together 
the stochastic nature of quantum systems (derived from the projection postu- 
late) with the usual unitarity evolution. In this process it is straightforward to 
convince oneself that the original Dirac's argument (cf. 0) is still valid, i.e. 
any linear superposition of states must be conserved during the evolution of the 
system, hence a linear operator can be defined which carries the initial states 
onto the evolved states at time t. 

In QED it is the virtual quantum fluctuations which drive the system from 
initial excited states into its ground state. This process is equivalent to rec- 
ognizing a patent stochasticity in the whole system (what other meaning can 
fluctuations bear?), a system which is necessarily compound (atom+field), even 
though there's no photon present in the field. Is it not much easier and more 
natural to assume that quantum systems possess an inherent random nature in- 
dependent of outer entities (fields, environments, etc)? Instead of rooting from 
a fundamental point of view the stochasticity of quantum systems in external in- 
controllable perturbations (observations, measurements), cannot we assert that 
this randomness is due to their quantum nature itself? Another way of stating 
this hypothesis is by realizing that in orthodox Quantum Mechanics, indetermin- 
ism is only present in the act of measuring, thus neglecting this indeterminism in 
the nature of the quantum systems itself (cf. page 108). As a matter of fact 
quantun systems which are not observed are deterministic. Nevertheless our 
experience in the laboratory clearly demands an indeterministic behaviour of 
these systems. Is it not more natural just to claim that these quantum systems 
are indeterministic themselves, with no resorting to external observations or in- 
terventions? Note that this does not suppose a major departure from quantum 
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orthodoxy, so we do not lose any result characteristic of Quantum Mechanics. 
Indeed the stochasticity, present in the formalism through the operator v(t) (cf. 
above), can be chosen to slightly modify the usual hamiltonian evolution. Sup- 
pose that the effects of the stochasticity are small compared to the hamiltonian 
evolution, i.e. we may write 



v' (s)v(s)ds 



1/2 



2 Jo 



v'(s)v(s)ds, 



(5.49) 



where v'(t)v(t) is of higher order than H in the norm sense. Then we may 
write for the density operator 



p(t)=PQM(t) + p st {t). (5.50) 

The usual well-known quantum results are contained in the first part, whereas 
new effects are only present in the stochastic part. As we have seen for two- 
level quantum systems, these new effects amount to introducing either energy- 
preserving terms which only show decoherence (only achieved by resorting to 
other systems in orthodox Quantum Mechanics) or energy- nonpreserving terms 
which in a natural way show spontaneous decay phenomena also without re- 
sorting to external interventions. 



On the other hand a mathematical scheme accounting for the stochastic 
nature of quantum systems is provided, as we have implicitly done, by letting 



theorem 2.1 be applied to each evolution operator component. Several com- 
ments should be remarked. First an elaborate argument should be offered for 
the cited theorem to be applied. This argument reads as follows. The stochas- 
ticity of quantum systems, as it is expressed by the projection postulate, roots 
on the impossibility of determining the evolution of a quantum system at a 
given instant t (typically when a measurement upon it is performed). At most 
we can know the probabilities of evolution from that state to the subsequent 
ones (typically the eigenvectors of the observable to be measured). The need 
for the presence of an observer has been the origin of neverending discussions 
since its very introduction. The mathematical way to deal with this kind of 
stochastic objects in general is to use random variables, which are nothing more 
than deterministic (usual) variables with a probability measure associated to 
them. So cannot quantum states, i.e. Hilbert-space vectors be also attached 
probability measures? That's all we have done. Once we have assumed such a 
hypothesis, the foregoing construction rests on mathematical results. If a linear 
operator tranforms quantum states into quantum states, then a stochastic linear 
operator will also transform stochastic quantum states into stochastic quantum 
states, hence endowing operators with an intrinsic random linear nature. This 
randomness can be expressed, as we have done, componentwise, so turning each 
operator (complex) entry into a complex random variable. Next step was then 
to apply theorem 2.1 to these complex random variables. 
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But a subtlety must be considered. When the whole process is built com- 



ponentwise (cf. eq.(3.30)), one may (and should) inquiry which basis is to be 
used? This surpringly drives us to the preferred basis problem. But now instead 
of trying to provide more or less complicate theoretical arguments, we may al- 
low Nature to freely choose the basis. To know what basis Nature chooses we 



make use of Prop. 4.1 to realize that the only basis (up to basis changes) in 
which the spontaneous decay is possible is the one which diagonalizes the free 
Hamiltonian. This observation follows from noticing that as it is usually done in 
a phcnomenological, though correct spirit to obtain decay we must add damp- 
ing terms — 7p. x to the time evolution equation of each of the density matrix 
components p. x expressed in the free hamiltonian eigenvector basis (cf. JI7J, 
page 64), which is equivalent to using a Lindblad operator of the form 
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Following Prop. 4.1 this Lindblad operator can only be obtained when using 
the previous stochastic formalism if H is diagonal, that is, in the free hamilto- 
nian eigenvector basis. Another way of stating it is by noticing that it is the 
spontaneous emission which chooses the basis. 



Finally the same hypotheses assumed in the orthodox quantum formal- 
ism (trace conservation, probability interpretation) are also taken up when the 
stochastic nature of quantum systems is considered. 

Note that the whole scheme is only based on a simple hypothesis: inherent 
randomness, the quantitative consequences of which are derived using just well- 
established physics-free stochastic methods. 



6 Conclusions 

In this paper we have considered the hypothesis of letting quantum systems 
show an inherent random nature. Mathematically we express this fact using a 
well-established theorem of decomposition of random variables using stochastic 
calculus. We then apply the usual definition of density operator for quantum 
systems being represented by stochastic state- vectors and impose normalization. 
As a result of these hypotheses we show how this stochastic evolution is deter- 
mined by the usual hamiltonian H of the system and a set of new operators 
Wfe(t) which characterizes the random behaviour of the evolution. The effect 
of these new elements in the description of the evolution of a quantum system 
makes possible to connect several originally independent facts. 

Firstly under suitable physical conditions, we may establish the general form 
of a completely positive Markovian evolution, i.e. a Lindblad evolution, the 
Lindblad operators being a more or less complicated expression of the hamilto- 
nian H and the operators Vk(t). Secondly, choosing adequately the operators 
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Vk (t) we may reproduce the phenomenon of spontaneous emission for a two-level 
quantum system in QED vacuum. Thirdly, choosing again these operators in a 
convenient way a system can show intrinsic decoherence, thus reducing the role 
of the observer in the measurement process. Finally, the formalism provides 
a natural way to attack the preferred basis problem, since now we have new 
phenomena such as spontaneous emission entering the physical picture which 
suggests a basis singularization. 

Nevertheless some remaining points should also be addressed. First and 
most important the question of a physical principle which provides the form 
of the operators Vk(t) must be investigated. In particular, in the same way as 
symmetry considerations usually provide a hint of the form of the hamiltonian 
H of a system, it would be desirable to have some sort of similar reasoning 
to work out the form of Ufe(t) in each concrete situation. Second, though the 
principles of the scheme are well established above, a generalization to more 
than 2 energy-level systems is necessary. Work in both these directions is in 
progress. 
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A Proof of Prop. g7j] 

The proof of proposition |2.l| is a conjugation of Ito's formula and the polar 



decomposition of operators acting upon Hilbert spaces. Using (|2.2j ) in (2.4) and 
making use of Ito's formula we get 

P (t) = EU st (t) P (0)EUl t (t) + 

+ a nm \ v n (s) p(0)v} n (s)ds , (A.51) 

i JO 

nm=l 

where the correlation matrix EfW^Wj™*! = a nm t has been taken into account. 



Now forcing trace preservation tvp(t) = trp(0) and using the cyclic property 
of the trace we may write 



ti[EUl t (t)EU st (t)p{0) + 

N 2 ft 

^ a nm / vl n (s)v n (s)dsp(0)] = tr 



■p(o), 

(A.52) 
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An = -K21I 2 


A12 


= ih\ 2 + 


|[(24i- l 22 )t 2X -l 


12^1] 


Ai 3 = -ih 2 i + 


^(2^1-^22)^12-^11] 


A14 


= K21I 2 




A 2i = 


--ih 21 + \{h 1 t 12 -l 21 t 22 ) 


A 22 


= -i{h\\ 


- h 22 ) + i[2£ n £* 22 - 


" Si=l 




A 23 = l 2X t 12 


A 24 


= -ih 2 i - 






A31 = 


^12 + 5^12-^22] 


A 32 


= ^12 






33 = i(hn - h 22 


) + i[2^ 22 -Eti l^l 2 


A34 


= ihi 2 + 


|Hti4 2 +i* 21 i 22 ] 






a 4 i = K12I 2 


A 42 


= -ihxi - 


f ^[4 2 (2^2-^i)- 


^2/22 


A43 = ih 2 i + 


WU^22-kl)-hlt 22 ] 


A44 


= -|4 2 | 2 







Table 3: Matrix entries A^'s in L[E n ] = X) m=1 A nm E m . 
which if it is to be valid for all p(Q) implies 

EUl t (t)EU st (t) + V a„ m / «t,(s)» n (s)ds = /. (A.53) 

nm=l ^ 

From this and making use of the polar decomposition of E[/ St (t) , we conclude 

/ * ft \ 1/2 

EU st (t)=U(t) \ i- J! 

(s)ds , (A.54) 

\ nm=l J ° J 

where U (t) is an arbitrary unitary operator which is expressed in the usual 
exponential form. 

B Proof of Prop. H3 

The idea of the proof is to use the linearity of the space A4 2 (C). We will 
firstly focus upon the canonical basis {En, E\ 2 , E 2 %, E 22 } which we shall rename 
respectively as {Ei,E 2 ,Ez,Ei\. Then we calculate L[Ek] for each fc. Denoting 

4 

L[E n ] = S A nm E m , 

m—l 

the results are shown in table ||. 

We will agree that these 16 quantities configure the 4x4 matrix A. By 
virtue of the linearity of M 2 (C) , we write 
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X = x n E n L[X) = x n L[E n ), 



(B.55) 



n=l 



n=l 



Then the condition L[X] = transforms into 



Oi x 2 x 3 x 4 ) ■ A= (0 0). 



(B.56) 



Now we must restrict the form of X, since it must be of the type X = vpov\ 
where po is a density matrix (positive self adjoint matrix with unity trace). Using 
again the linearity of M. 2(C), we may express 



4 

n=l 



Then X may be written as X — ^ n a n E n with 



Oil 



= (Px fo)p°(f*) "2 = (Pi &)a>(^' 



"3 = (p3 At) PO ( ^* J "4 = (03 Pi) PO [pi 



These a^s show the following immediate properties 

1. ai = Vpo <^ /3i = Pi = => a 2 = Vp 

2. a 4 = Vp & Pz = Pi = a 3 = Vp 

3. Q!2 = V/9o «3 = Vpo 

4. a 2 = Vpo => /?i = P2 = or fi 3 = /3 4 = 0. In the first case we get 
ail = whereas in the second we get 0:4 = 0. 

Then, in matricial form, X = vpoV' may show the following aspects: 



• liv 







, then X = ol\E\. 



• If v = I ? ? 1 , then X = 014E4. 

k P3 P4/ 

• If w = ( ^ ^ 2 I with at least three /3t's not equal to 0, then X 



En=l a n,E n 



Carrying these a^s to (B.56) and taking into account the expression for A^ 's 
we obtain the conditions announced in Prop. 4.1. 
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C The Jaynes-Cummings model 



The Jaynes-Cummings hamiltonian can be straightforwardly diagonalized by 
noticing that in the |±,n) basis it adopts the matricial representation 



/ - wo/2 
' uj /2 





c 

-loq/2+uj 





oJo/2+noj 
fVn+1 -wo/2+(n+l)o; 



V 



II 



(I 



from which we extract the eigenvalues (n > 0) 



E- 



(n+-)w± 
-wo/2 



W — CiJQ 



+ (n+ l)e 2 



and the eigenvectors (n > 0) 



K+i> 

\ u o) 



cos n+1 \ 
- sin 6 n +i\ 

-,o), 



,n) + sin0„ + i|-,n + 1) 
, n) + cos0„+i| — , ri + l) 



where 



tan 6>, 



2 



n+1 



+ ^/( = ) 2 + (n+ l) e 2 



£\/n + 1 



These are the eigenvalues and eigenvectors of the dressed atom, which is 
commonly used in quantum optics (cf. fujfl ). Notice that the set \uq )} 

constitutes an orthonormal basis for the joint Hilbert space of atom+field. 
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